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La pr6sente invention conceme une m6thode pour compacter et filtrer des donn6es 
g6ophysiques dans le but d'en extraire au mieux des informations sur la nature du sous-sol. 

Cette raethodologie s'applique i la compaction et au filtrage des jeux de donnees 
10 geophysiques multidomaines, i.e. soit des donn6es sismiques avant sommation issues 
d'une campagne sismique 3D regroup6es en un certain nombre de cubes de donn6es • 
correspondant chacun par exemple h un mSme deport (cubes de traces sismiques dits iso- . 
offsets) ou a un meme angle d'incidence (dits cubes iso-angles dMncidence). ou h des cubes 
d'attributs sismiques obtenus par traitement des amplitudes sismiques comme. par .< 
15 exemple. des cubes de propri6t6s elastiques issus d'une inversion stratigraphique jointe de 
donndes sismiques avant sommation, ces cubes 6tant obtenus par la mise en oeuvre de la" 
m6thode d^crite par exemple dans le brevet FR 2 800 473 (EP 1 096 270) du demandeur. H 
s'agit alors de mettre en forme rinfonnation sismique avant sommation contenue dans des 
portions de traces sismiques formes par exploration d'une zone souterraine, k I'aide de 
nouveaux attributs repr6sentatifs qui sont ensuite exploit6s pour I'inte^pretation geologique 
ou rinterpr6tation dhm reservoir du sous-sol. Hie peut aussi s'appUquer & la compaction et 
au filtrage de donnees sismiques 3D avant ou apr^ sommation de campagnes r6p6t6es 
dans le temps (sismique 4D), ceci afm d'interpr6ter les variations de contenu en fluides 
d'un reservoir p6trolier en cours de production. Cette m6thode peut aussi s'appUquer a la 
25 compaction de donn6es diagraphiques de plusieurs natures physiques, en formant k partir 
des mesures de d6part de nouveaux attributs diagraphiques, et ce afin d'am^liorer 
I'inteipretation des donn6es diagraphiques en termes de caracteristiques lithologiques et 
p6trophysiques du reservoir. 
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Etat de la technique 

Les mesures geophysiques sont classiqueraent utilisees pour foiimir des informations 
<:ompl6mentaires, par rapport aux observations directes dans les puits, sur les variations des 
proprietes geologiques des fomaations du sous-soL L*interpr6tation sismique, par exemple, 
est bas6e sur 1* analyse d'attributs aiSn d'extraire ces propri6t6s g6ologiques & partir de 
r information sismique, et ce dans une zone d'interSt d'une formation souterraine (zone 
reservoir par exemple). Un autre exemple est Tinterpretation diagraphique oii les diverses 
mesures sont analysees simultanement pour en d6duire les caracteristiques petrophysiques 
de la formation souterraine etudiee. 

Une approche standard pour T analyse d'attributs sismiques ou diagraphiques 
implique T utilisation de m^thodes statistiques de reconnaissance de forme et d'estimation 
afin de calibrer les attributs pertinents avec Tinfoimation disponible obtenue au droit de 
differents puits fores. 

Par exemple, par le brevet FR -A 2 768 818 (US 6 051 651 du demandeur, on connait 
une m6thode permettant de detecter les classes potentielles dans une population 
d'evenements li& aux proprietes physiques d'un milieu complexe tel que le sous-sol, 
rep6res k partir de donnees obtenues par exploration du sous-sol ou par des mesures in situ, 
(ev6nements d'ordre ggologique, g6ophysique, evenements lids k des donnees de 
production de fluides etc.), ces 6v6nements etant reperes par des points dans un espace 
multivariable et defmis par un vecteur d'attributs. 

Un autre exemple est donne par le brevet FR 01/05.675 du demandeur, qui d6crit une 
mdthode pour faciliter le suivi au cours du temps de revolution d'une zone souterraine par 
analyse comparee d*un certain nombre n de jeux d'enregistrements sismiques obtenus 
successivement k Tissue de n campagnes sismiques 3D successives (sismique dite 4D), 
dans laquelle on utilise une technique de reconnaissance des foraies appliqu6e h Tensemble 
des evenements sismiques de plusieurs campagnes, considdres et analys6s simultanement. 

Citons aussi le brevet EP -A-671 017 (US 5 638 269) du demandeur, qui d6crit une 
methode permettant de rapprocher des donnees g6ologiques obtenues par carottage ou 
diagraphies dans des puits et des donnees sismuques obtenues par une campagne 
d'exploration sismique qui s'appuie sur une technique de calibrage statistique avec un 
6talonnage statistique obtenu en rapprochant les donnees geologiques locales (mesurees 



dans des puits) avec des attiibuts sismiques lus sur des traces sismiques obtenues au 
voisinage imm6diat de chaque puits. 

Les attributs sismiques qui font I'objet de ce type d'interpretation sont calcules dans 
le domaine apres sommation, ou dans le domaine avant sommation. 

5 Dans le domaine apits sommation, les attributs classiques sont calcul6s ^ partir des 

amplitudes au niveau du reservoir, ou k partir de I'imp^dance P estimee k I'aide d'une 
inversion stratigraphique. 

Dans le domaine avant sommation, le nombre d'attributs pertinents peut 
considerablement augmenter. Dans ce cas, les attributs sont soit les amplitudes k diffdrents 
10 offsets (deports) ou angles d'incidence, soit des parametres issus d'une inversion 
stratigraphique jointe avant sommation (Impedances P et S, masse specifique, produit des 
parametres de Lam6 avec la masse specifique, etc.). 

En ce qui conceme I'interpretation des mesures diagraphiques, la situation est 
analogue, un trds grand nombre de mesures etant disponibles d'ovl en decoulent un grand 
15 nombre d'attributs permettant de caracteriser les diagraphies, ces attributs 6tant souvent 
organises en families h mSme signification physique (par exemple, famille des attributs de 
r6sistivit6, famille des attributs associds k la radioactivity natureUe ou induite des 
formations, etc.). 

Une mdthode dite d'analyse en composantes principales (ACP) bien connue des gens 
20 de I'art, peut Stre appliqu6e k ces attributs afin d' analyser les relations existantes entre eux 
et de reduire le nombie d'attributs. significatifs. Les composantes principales extraites de 
I'ACP sont les nouveaux attributs : ils d6fmissent une base vectorielle orthogonale ou non 
(rotation) et sont des combinaisons lin6aires de variables initiales. lis peuvent etre utilisees 
dans des algorithmes statistiques de reconnaissance de forme et correspondent k un filtrage 
25 multivari6 de I'information sismique ou diagraphique initiale. Diff6rentes applications de 
I'ACP sont d^crites, par exemple, dans : 

- Dumay, J., Foumier, F., 1988, "Multivariate statistical analyses applied to seismic 
facies recognition", Geophysics, S3, 1151-1159 ; 

- Hagen, D., C, 1982, "The application of principal components analysis to seismic data 
30 sets" : Geoexpl., 20, 93-111 ; ou encore 




- dans le brevet FR-A- 2 772 138 (US 6 345 108) du demandeur. 

Mais TACP ne prend pas en compte la notion de groupe d'attributs et, par 
consequent, il est souvent difficile de donner une interpretation physique aux composantes 
principales, ou bien de les relier clairement aux attributs initiaux, surtout si ceux-ci sont 
5 organis6s en groupes de meme signification physique, 

L' analyse dite canonique est aussi une methode statistique bien connue des gens de 
Tart, qui permet d'6tablir les relations pouvant exister entre deux groupes de variables afin 
de savoir si ces deux groupes decrivent les memes propri6tes. Cette methode est decrite par 
exemple par : 

10 - Hotelling, H., 1936, "Relations between two sets of variables", Biometrika, 28, 321- 
377. 

Un exemple d* application de 1' analyse canonique est d6crite par exemple par : 

- Foumier, P., and Derain, J.F., 1995, "A statistical methodology for deriving reservoir 
properties from seismic data". Geophysics, 60, 1437-1450. 

15 Cette m6thode se limite a T etude de deux groupes de variables et done ne s' applique 

pas de fagon g6n6rale k Tinfomiation sismique multicubes ou k Tinformation diagraphique 
multidomaines. De plus, elle d^finit des variables synlii6tiques dans chacun des sous- 
espaces associ6s aux deux groupes, et non une seule base vectorielle permettant de 
compacter 1' ensemble des variables initiales et de d6crire les relations entre les deux 

20 groupes. 

Plusieurs m6thodes permettant de g6n6raliser Tanalyse dite canonique en 6tendant 
r analyse a plus de deux groupes de variables existent. DiffSrents aspects de Tart anterieur 
dans le domaine consider^ sont decrits, par exemple, dans les publications suivantes : 

- Horst, P., 1961, "Relations among m sets of measures", Psychometrika, 26, n*2, 129- 
25 149 ; 

- Carol!, J.D., 1968, "A generalisation of canonical correlation analysis to three or more 
sets of variables", Proc. 76* Conv. Amer. Psych. Ass ; 

Kettenring, J.R., 1971, "Canonical analysis of several sets of variables". 



58, 3, 433-450 ; ou 

- Saporta, G.. 1990, "Probabilit6s. analyse des donnees et statistiques", Technip, Paris. 

Ces m6thodes pr6sentent certaines limites. En particuUer. elles ne pennettent pas de 
d6crire chaque groupe separement, les variables synth6tiques ne peuvent 6tre reU6es 
facilement aux diffdrents groupes. et . par cons6quent, 6tre physiquement inteipt6t6es. De 
plus, seule la part d'inertie globale qu'elles represented pent 8tre connue, et non la part 
d'inertie qu'elles lepr^sentent pour chaque groupe. D n'est par ccns6quent pas possible de 
filtrer chaque groupe ind6pendamment. 

Une technique d'analyse dgalement connue dite analyse en composantes principales 
g6n6ralis6e (ACPG). peimet de comparer diff6rents groupes de variables tout en d6crivant 
chacun d'eux. On la trouve mise en ceuvre sur des donnees i caractfere economique par 
exemple dans : 

- Casin, Ph. 2001, "A generalization of principal component analysis to K sets of 
variables". Computational Statistics & Data Analysis, 35, 417-428. " 



La m6thode sdon rinvention 

La m6thode selon I'invention propose une approche ot en exploitant une technique 
statistique multivariable appliqude k plusieuis jeux d'enregisttements sismiques ou 
diagraphiques disponibles. obtenus par exploration tfune zone soutcrraine, on compacte et 
on filtre les ev6nements g6ophysiques lisibles sur ces jeux de mesures. La m6thode 
comporte : 

- la repartition de ces ev6nements g6ophysiques dans un certain nombre p de families 
d'attributs g6ophysiques ayantchacune une signification physique particuli^re ; 

- la formation par combinaison des dits attributs g6ophysiques, de variables synth6tiques 
en nombre trfes inf6rieur au nombre des dits attributs. que Ton obtient par construction 
d'une base vectorielle orthogonale dans chacun des p ensembles d'analyse constitu6s 
par les donn6es de chacune des dites families. d'oCi d6coule la formation d'une base 
vectorielle orthonorm6e d6crivant les p ensembles d'analyse ; 
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- I'utilisation de cette base vectorielle orthonormee pour filtrer et deciire les dits 
6venenients g6ophysiques. 

La formation des variables synth6tiques est obtenue par exemple en combinant une 
technique d'analyse en composantes principales ou ses techniques d6rivees bastes sur la 
diagonalisation d'une norme (recherche de vecteurs propres) et une technique de 
regression simple. 

Suivant un mode de mise en oeuvre, on peut soxmniettre les variables synth6tiques a 
une rotation pour ameliorer les correlations entre variables initiales, c'est a dire les attributs 
geophysiques, et variables synth6tiques, 

Les variables synth6tiques peuvent etre form6es k partir de families comportant par 
exemple un ou plusieurs cubes de traces sismiques obtenues k Tissue d'une campagne 
sismique 3D, chacun correspondant k un meme deport ou offset. 

On peut ^galement former des variables synth6tiques k partir de families constituees 
de cubes de traces sismdques obtenus k Tissue d*une campagne sismique 3D, chacun 
correspondant h un meme angle d*incidence. 

On peut encore former des variables synthetiques k partir de families constituees de 
cubes d* attributs sismiques formes par traitement des amplitudes de traces sismiques 
obtenues k Tissue d'une campagne sismique 3D. 

On peut encore former des variables synthetiques k partir de families constituees de 
cubes de traces sismiques obtenues k Tissue de campagnes sismiques successives de la 
zone. 

On peut encore former des variables synthetiques k partir de families constitutes de 
cubes d* attributs sismiques form6s par traitement des amplitudes de traces sismiques 
obtenues k Tissue de campagnes sismiques successives de la zone. 

On peut encore former des variables synthetiques k partir de families constituees de 
mesures d'attributs diagraphiques de nature physique differente, comme la famille des 
mesures de r6sistivit6, la famille des mesures li6es k la radioactivite des terrains, la famille 
des mesures li6es k la porosit6 des terrains, etc. 
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La methode est particulierement avantageuse en ce qu'elle permet, appliquee k de 
• rinformation sisinique ou diagraphique: 

- de definir de nouveaux attributs geophysiques, combinaisons lineaires des attributs 
gdophysiques initiaux, r6sumant rinfonnation sismique ou diagraphique initiale et 

5 permettant de filtrer cette information, tout en respectant la structure multidomaines 

des donnees 

- de montrer si toutes les families d' attributs geophysiques initiaux, ou seulement 
certaines d'entre elles, apportent une information redondante, et de mesurer la part de 
cette information redondante 

10 - d' apporter une interpretation physique aux nouveaux attributs geophysiques. 

Les nouveaux attributs g6ophysiques ainsi formes peuvent ensuite etre analys6s afin 
dextraire les proprietes g6ologiques de Hnfonnation sismique ou diagraphique. Deux 
classes de m6thodes peuvent classiquement etre utilis6es k cette fin : une premiere qui 
utilise des techniques d'estimation statistique et foumit une interpretation quantitative en, 
15 termes de propri6tes geologiques, et une deuxi^me qui utilise des techniques de 
reconnaissance de forme pour analyser les facifes sismiques ou les facies diagraphiques 
(electrofacids) et foumir une interpretation qualitatitive. 

Des exemples de ces deux types de m6thode sont xms en oeuvre respectivement dans 
les brevets EP~A-671 017 (US 5 638 269) et FR-A-2 768 818 (US 6051651) du 
20 demandeur, d6j& cites. 

Presentation succincte des figures 

- les figures lA a IC montrent 1' organisation des enregistrements sismiques en p 
groupes : chaque groupe rassemble dans un cube T ensemble des portions de traces 

25 correspondant k un critfere physique identique ; 

- les figui-es 2A k 2C montrent respectivement les cartes de trois premiers nouveaux 
attributs sismiques Z^^\ Z^^^ et Z^^^ ; 

- les figures 3A k 3D montrent respectivement le premier plan temporel vrai (3A) de la 





fenetre d'analyse et le m6me plan reconstitu6/filtre avec 1 (3B), 3 (3C)ou 8 (3D) des 
nouveaux attributs sismiques pour le cube d' impedance P ; 

les figures 4A a 4D montrent respectivement le premier plan temporel vrai (4A) de la 
fenetre d'analyse et le mame plan reconstitu6/filtr6 avec 1 (4B), 3 (4C)ou 8 (4D) des 
nouveaux attributs sismiques pour le cube d'impedance S ; 

les figures 5A a 5D montrent respectivement le premier plan temporel vrai (5A) et le 
m6me plan reconstitu6/fOtr6 avec 1 (5B), 3 (5C)ou 8 (5D) des nouveaux attributs 
sismiques pour le cube de miasse specifique p ; 

les figures 6 montrent respectivement les sections verticales. passant par le puits W2, 
vraie (6A), reconstitu6e/filtr6e avec les 8 premiers nouveaux attributs (6B), ainsi que 
les differences (i^sidus) (6C) entre ces deux sections pour le cube d'imp6dance P ; 

-les figures 7 montrent respectivement les sections verticales, passant par le puits W2, 
vraie (7 A), reconstim6e/filtr6e avec les 8 premiers nouveaux attributs (7B), ainsi que 
les differences (r6sidus) (7C) entre ces deux sections pour le cube d'imp6dance S ; 

la figure 8 montre la localisation des traces d'apprentissage utilis6es pour calibrer la 
fonction de classement dans le cadre de la reconnaissance des formes supervisee des 8 
premiers nouveaux attributs sismiques ; 

les figures 9A. 9B montrent respectivement les attributions des ev^nements sismiques 
et les attributions les plus fiables de ces mgmes 6v6nements dans le cadre de la 
reconnaissance des foraies supervisee des 8 premiers nouveaux attributs sismiques ; et 

• les figures lOA, lOB repr6sentent respectivement les attributions des 6venements 
sismiques et les attributions de ces mSmes 6v6nements sismiques dans le cadre de la 
reconnaissance des formes non supervisee des 8 premiers nouveaux attributs sismiques. 



Description d€taillee 

On dispose de plusieurs jeux d'enregistrements g6ophysiques obtenus par exploration 
sismique ou reconnaissance diagraphique d'une zone souterraine 6tudi6e. Les 6v6nements 
geophysiques que Ton cherche \ compacter/filtrer sont caract6ris6s par des attributs 



geophysiques souvent organises en groupes ou families a signification physique 
particuliere (figures lA k IC). Ces attiibuts peuvent 6tre la succession des valeuis des 
amplitudes de portions de traces sismiques ou de paramdtres 61astiques issus d'une 
inversion elastique jointe avant sommation, obtenue par exemple & I'aide de la mfithode 
d6crite dans le brevet ER 2 800 473 (EP 1 096 270)d6ja citd le long de la fen6tre 
temporelle 6tudi6e. Ces attributs peuvent gtre, pour une profondeor donn6e, les mesures 
d'un mSme ensemble d'enregistrements diagraphiques realis6s le long d'un puits, comme 
par exemple les valeurs des diff6rents outils de r6sistivit6, les valeurs des diff^rents outils 
mesurant la porositd, etc. 

1 - Notations 

Soit X une matrice compos6e de p tableaux juxtaposes Xj, , . . Xj . . . Xp : 
X ={X. ... X, ... X,} , Xj=[X,, ... X,, ... X,„^] 

oii Xj est un tableau compose de mj variables Xj.k de n individus (au sens statistique-?du 
terme). La matrice X rassemble tous les 6v6nements g6ophysiques disponibles. Ces 
dv^nements sont tout d'abord organises en p groupes ou families Xj (figure lA). Chaque 
faraille Xj correspond k une propri6t6 physique particulifer©. Pour I'analyse de la sismique 
3D, cette famille se pr6sente sous la forme d'un cube de valeurs correspondant k cette 
propri6t6, rep6r6es par leur position spatiale (ligne, cdp, temps/profondeur), (figures IB et 
IC). Ces cubes peuvent Stre formes par exemple par un ensemble des valeurs d'amplitude 
de n traces ou portions de traces associ6es k un d6port conmiun ou k un angle commun ou 
encore la distribution spatiale des valeurs d'un attribut commun d6riv6 de ces amplitudes 
telle que, par exemple, l'imp6dance acoustique en ondes P ou S. la masse specifique ou des 
paramfetres de Lam6, dans une fenStre temporelle choisie de taille mj. Dans cet exemple, 
chaque variable Xj,k rassemble les valeurs de I'attribut sismique consid^re (amplitude 
sismique ou propriet6 Elastique) pour un echantillon temporel k Tintdrieur de la fenetre 
d'analyse. Dans cet exemple, un cube comprend autant d'individus que de traces 
disponibles. Pour I'analyse des enregistnements diagraphiques, cette famille se pr6sente 
sous la forme d'une suite de mj s6ries de valeurs, repdi^es par leur position le long du puits 
(temps ou profondeur). Cette suite peut 8tie form^e, par exemple. par I'ensemble des 
mesures diagraphiques des outils de mSme nature physique, comme les diff6rents outils 
mesurant la porosit6 de la foraiation, ou les outils mesurant diff^rents types de r6sistivite, 
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sur un ensemble de n cotes pxofondeur ou temps le long d'un ou de plusieurs pmts" au 
niveau de la formation souterraine 6tudi6e. Dans cet exemple, chaque vanable X,, 
rassemble les valeurs de Tattribut diagraphique consider^ (r6sistivite, porosite, 
radioactivit6. ...) pour m.e cote profondeur ou temps ^ rint6rieur de la fenStre d'analyse. 
Dans cet exemple. une suite comprend autant dMndividus que de cotes profondeur ou 
temps disponibles le long du ou des puits analyses. 

Les m, variables Xj.. d'une table Xj definissent un espace d6not6 et sont supposes 
centr6es. Soit D„ la matrice diagonale des poids des individus. 

Soit Uj un op6rateur de projection orthogonale sur r espace Sj d6fmi par : 

2 - Calcul des nouveaux attributs 

Le but de r Analyse en Composantes Principales Gen6ralisee (ACPG) est I la fois de 
d^rire les proximit6s entre les p families de variables et de d6crire separement chaque 
groupe grtce a des variables synth6tiques qui sont des combinaisons Uneaires des vanables 
15 initiales. Alors, ce qui est commun 2l tous les groupes ou families, ou k certains d'entre eux 
seulement, et par cons6quent ce qui est diffdrent, pent 6tre mis en Evidence. L'ACPG 
pennet de r6duire le nombre de variables utilis6es pour dficrire un jeu de donn6es et permet 
un filti^ge multivari6 tout en conservant la structure multidomaines des donnees miUales. 
Afin d'atteindre ce but, la premiere variable synth6tique calcul6e, que I'on nomme Z , 
20 calcul6e doit Stre telle que : 

1. les p projections orthogonales Z,<» de Z<»> sur les espaces Sj doivent etre le plus 
proche les unes des autres ; 

2. Dans chaque groupe j la dispersion des n individus projetes sur z/^> est 
maximum. 
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1. ±R\Z'\Z^') 

2. imx |;Cov^(Z«^,J 

v=l 

Ceci se traduit mathematiquement par les conditions suivantes : 

ou est le coefficient de correlation multiple et Zj^^^ est la projection de Z^^^ sur Tespace 
Sj normee k 1. 



5 Ensuite, on cherche une nouvelle variable synthetique Z^^ selon les mSmes criteres. 

On ajoute une condition d'orthogonalite entre les projections les variables Z^^ du mSme 
groupe. Le processus s'axrSte quand une base orthogonale dans chaque espace Sj est 
trouvee. Le nombre de variables Z^^^ n'excede pas la dimension la plus grande des. espaces 

Sj. 

10 Des methodes de rotation bien connues des gens de Tart, peuvent 6tre mises en oeuvre 

dans chacun des espaces ainsi trouv^s, afin d*optimiser les coefficients de correlations 
entre les nouvelles variables et les variables initiales. 

La premifere variable Z5*^ est vecteur propre nonn6 de la matrice de covariance 
XX'Dn correspondant & la plus grande valeur propre. Z^^^ est ensuite projete sur les p 

15 espaces Sj et, dans chaque espace, la regression des mj vaiiables par Z/^^ est calcul6e. La 
matrice des p r6sidus juxtaposes Xj^^^ definit le nouvel espace d'6tude. Ce nouvel espace 
correspond ^ Tespace initial duquel Tinfluence de Zj^^^ a 6t6 retire dans chaque groupe Sj. 
Les raj variables de la table Xj^^^ definissent un sous espace Sj^^^ de Sj. La seconde 
composante Z^^^ est le vecteur propre norme de la matrice de covariance X^^^X^^^^n des 

20 residus correspondant k la valeur propre la plus grande. Z^^^ est projetee sur les p espaces 
Sf^\ Et la regression des mj variables par Zj^^ est calcul6e dans chaque espace Sf^\ La 
regression effectuee dans chaque espace Sj assure que la projection de Z^^^ sur Tespace Sj^^^ 
est orthogonale k la projection de Z^^^ sur Tespace Sj. 
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Le calcul des vecteurs propres est effectud par toute methode connue des gens de 



Tart. 



En pratique, la matrice XX'Dn 6tant d'ordre n. n pouvant gtre tr^s grand, k une 6tape 
k. on calcule le vecteur propre norme V^> de la matrice X<^^^D„X<^^ associe a la plus grande 
> valeur propre. Le vecteur est ensuite simplement calculi par transformation Uneaiie 
2;W=x<^)V*> et norme. 

Chaque groupe peut 6tie reconstitu6 avec un petit nombre d'attributs. Les variables 
initiales du groupe j peuvent s'exprimer comme la somme d'une combinaison Un6aire des 
variables 7^^'> et d'un r^sidu. Si toutes les variables Zj<^> d'un groupe sent utiUs6es, le 
0 r&idu est nul. Les coefficients de la combinaison lineaire ne sont autres que les 
coefficients calculus k I'issu de chaque regression dans les espaces Sj. 

La m6thodc selon invention a pour propriet6, entre autres, que les variables Z^^^ sont 
orthogonales devix k deux. 

En vue de I'interpretation des variables synthetiques calcul6es. on calcule ensuite 
15 . pour chaque groupe j : 

1. le carre de la corr61ation entre la variable Z*^ et la variable 7^^ du groupe j. 
indiquant le degre de proximite entre ces deux variables; 

2. le pourcentage d'inertie du groupe j expliqu6 par la variable Zj*\ 

Les correlations entre les variables Z*^ et les variables initiales X,, sont 6galement 
20 calculus pour mettre en Evidence d'6ventuelles relations entre les Xj., . l^s correlations 
entre les variables Zj® et Xjj. d'un m6me groupe j, ou les correlations entre variables Zj»> 
pour un mSme k peuvent 6galement 6tre calculees. 

Les nouvelles variables synth6tiques ainsi cr€6es constituent les nouveaux attributs 
g6ophysiques recherch6s, 

25 3 - Exemple d'application 

La m6thode selon I'invention a 6t6 appliqu6e k trois cubes de propri6t6s elastiques 
(impedance P not6e Ip. imp6dance S notee Is et masse specifique not6e p) obtenus apr^s 
inversion 61astique jointe avant sommation k 1' aide de la methode decrite dans le brevet FR 




2 800 473 dejk cite d'enxegistrements sismiques regroupes en cubes k angles dMncidence 
commun. Chaque propriete elastique engendre un groupe de variables ou d'attributs 
sismiques. La fenetre temporelle d' analyse se compose de 19 echantillons de signaux 
pr61ev6s sur les enregistrements sismiques. Nous consid6rons comme variables les valeurs 
5 de la propri6t6 elastique pour un 6chantillon temporel a Tinterieur de la fenetre d' analyse, 
Chaque groupe est alors compos6 de 19 variables composees d'autant d^individus que de 
traces disponibles. Comme il y a trois groupes de variables, nous disposons au total de 57 
variables. Les trois cubes de paramdtres ont 6t6 pr6alablement horizontalises par rapport ^ 
une surface de r6ference. 

10 Les variables synthetiques Z^^^ obtenues r6sument les informations respectives des 

trois proprietes elastiques, tout en prenant en compte les relations entire elles. 

3.1 - Analyse des relations entre les impedances P, S et la masse specifique p 

Le tableau 1 ci-dessous donne la proportion de variance en pourcentage expliqu6e 
par les projections des 8 premieres variables synthetiques Z^^^ h Z^^^ dans les trois 
15 groupes. II montre quMl suffit de 8 variables pour resumer 87,4% de la variance de 
rimpedance P, 94,4% de la variance de 1* imp6dance S et 88,4% de la variance de la masse 
sp6cifique. 







2(2) 




Z(4) 


Z(5) 


2^(6) 


2P) 


2;(») 


• 


Ip 


23,4 


24,5 


6,3 


13,1 


9,8 


3,1 


3,8 


3,4 


87,4 


Is 


26,6 


14,2 


15,2 


12,6 


5,7 


7,2 


5,8 


7.1 


94,4 


• 


16,9 


14,1 


13.6 


10,2 


8,0 


9,4 


10,1 


6,1 


88,4 



20 Le tableau 2 ci-apr&s donne le carr6 de la corr61ation entre les variables synthetiques 

Z^^^ et leurs projections sur les trois groupes, Cette valeur indique le degr6 de proximit6 
entre la variable et le groupe. Plus la valeur est proche de 1, plus la variable est proche de 
sa projection et done du groupe. 
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2,(1) 


2^(2) 


2,(3) 


2;(4) 


2(i» 


2^(0) 




2(«> 


Ip 


0.48 


0,79 


0,38 


0,72 


0,71 


0.24 


0,02 


0,30 


Is 


0,91 


0,43 


0,86 


0.53 


0,36 


0.59 


0.09 


0.73 


• 


0.78 


0,49 


0,77 


0.46 


0,65 


0,79 


0.96 


0,46 



Par exemple. Z<'> est tifes proche de Is et proche de p , mais ties eloigne de Ip. Mors 
que Z«> est proche de Ip, mais trfes 61oigne de Is et de /? . Sur les 8 premiferes variables Z»\ 
trois sont proches de Ip seul, deux sont proches ^ la fois de Is et dep , une est proche de Is 
seal, et deux sont proches de p seul. Ainsi, 11 apparait que Ip et Is ne contiennent pas la 
m6me information et qu'une part de rinfoimation contenue dans Is et dans p semble 
conunvme. 

3.2 - Filtrage des impedances P, S et de la masse specifique 

Les trois propri6t6s 61astiques peuvent 6tre reconstitutes avec un nombre restreint de 
nouveaux attiibuts sismiques. Le nombre de variables synth6tiques utilis6es pour la 
reconstitution pent Stre diff6rent d'un groupe k I'autre : il depend du pourcentage d'inertie 
que Ton d6siie reconstituer/filtrer. 

3.3 - Reconnaissance des formes supervisee et non supervisee des nouveaux 
attributs 

II est possible d'interpreter geologiquement les nouveaux attributs sismiques au 
travers d'une analyse de fades avec un algorithme de reconnaissance de forme supervise 
ou non. 

En reconnaissance des formes supervisee, il est necessaire de choisir des classes 
d'apprentissage regroupant des 6v6nements sismiques situ6s aux voisinages de puits dont 
!■ interpretation g6o]ogique a et6 effectu6e. Cette m6thodologie permet de classer 
I'ensemble des 6v6nements sismiques en fonction des caracteres g6ologiques observes aux 




puits. 

Dans notre cas, les connaissances geologiques sur les trois puits disponibles Wl, W2 
et W3, ainsi que les connaissances geologiques generales sur la zone d'6tude nous amenent 
& considerer trois classes : une premiere classe correspondant aux puits Wl et W2 
5 representent des sables turbiditiques, une seconde classe correspondant au puits W3 
represente des boues turbiditiques et une troisieme classe repr^sente des hemipelagites. Les 
puits Wl et W2 montrant des caracteristiques geologiques proches, seuls les 6v6nements 
sisnniques situes autour du puits W2 sont utilises, le puits Wl servant de test aveugle. La 
localisation du puits W3 etant incertaine, les ev6nements sismiques utilises pour 

10 caracteriser la classe 2 sont situ6s autour d'un puits virtuel defini k proximit6 de W3. 
Enfin, ne disposant d'aucun puits pour la classe 3, un puits virtuel est d6fini k partir de nos 
connaissances geologiques g6n6rales. La Fig.7 montre la localisation g6ographique des 
trois classes d'apprentissage. Chaque classe est constitute de 121 ev6nements sismiques 
(11x11). Ces 6v6nements servent a calibrer une fonction de classement qui .fest ensuite 

15 appliquee k T ensemble des 6venements sismiques. 

Les Fig. 8 A, 8B montrent respectivement les affectations obtenues et les affectations 
les plus probables, c'est-2t-dire ayant une probabilite d'affectation superieure ou egale k 
80%. Le puits Wl, utilis6 en test aveugle, est coirectement attribu6 a la meme classe que le 
puits W2. Le puits W3 est localise en bordure des classes 1 et 2. La classe 2 se distribue le 
20 long d'une forme chenalisante. La carte de la figure 8B montre que les r6sultats sont 
stables : peu d*evenements presentent une mauvaise probabilite d' affectation. 

En reconnaissance des formes non supervisee, les classes d'apprentissage sont 
constitutes en utilisant les modes d'une fonction de densitt de probabilite multivariee 
calculee k partir de Tensemble des evtnements sismiques, dans un espace engendrt par les 

25 huit premiers nouveaux attributs. Les 6v6nements sismiques de plus forte densite de 
probability constituant les pics de cette fonction de densite sont selectionn6s pour realiser 
Tapprentissage de la fonction de classement. Cette fonction est ensuite appliqute a 
I'ensemble des 6venements sismiques consider6s pour les affecter k Tune des classes : les 
Fig.PA, 9B montrent les affectations obtenues, ainsi que les affectations les plus probables, 

30 c'est-k-dire ayant une probability d' affectation sup6rieure ou 6gale k 80%. Les cartes 
obtenues sont trfes proches des cartes des fig. 8A, 8B. La classe 3 de Tanalyse supervis6e 
correspond aux classes 3, 4 et 5 de I'analyse non supervisee. Les puits Wl et W2 sont 
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attribu6s k la mame classc et le puits W3 est toujours situe k la frontiere entre les classes 3, 
et surtout 1 et 2. Cette convergence des i€sultats valide l'inteipr6tation de Vapproche 
supervisee. 




REVENDICATIONS 



1) Methode pour compacter et filtrer des ev€nements g6ophysiques lisibles sur 
plusieurs jeux d'enregistrements g6ophysiquesdisponibles, obtenus par exploration d'une 
zone souterraine, caract6risee en ce qu'elle coniporte : 

5 - la repartition de ces 6v6nements gfiophysiques dans un certain nonabre p de families de 
variables g^ophysiques ayant chacune une signilBcation physique particuliere ; 

- la fonnation par combinaison des dites variables sismiques, de variables synth6tiques 
en nombre trts inf6rieur au nombre des dites variables, que Ton obtient par 
construction d'une base vectorielle orthogonale dans chacun des p ensembles d'analyse 

10 constitutes par les donn6es de chacune des dites families, d'oil decoule la foraiation 
d*une base vectorielle orthonomi6e dficrivant les p ensembles d* analyse ; et 

- Tutilisation de cette base vectorielle orthononn6e pour filtrer et d6crire les dits 
6v6nements g6ophysiques. 

2) M6thode selon la reyendication 1, caract6ris6 en ce que la formation des variables 
15 synth6tiques est obtenue en Combinant une technique d' analyse en composantes principales 

ou des m^thodes d6riv6es de diagonalisation d'une norme et une technique de regression 
simple. 

3) M6thode selon la revendication 1 et 2, caracteris6e en ce que Ton soumet h une 
rotation les variables synth6tiques pour amSliorer les correlations entre variables initiales et 

20 variables synthetiques. 

4) M6thode selon la revendication 1 ou 2 ou 3, caracterisee en ce que Von forme des 
variables synthetiques k partir de families constitutes de cubes de traces sismiques obtenus 
a Tissue d'une campagne sisndque 3D, chacun correspondant & un meme deport ou offset. 

5) Methode selon la revendication 1 ou 2 ou 3, caracterisee en ce que Ton fomie des 
25 variables synthetiques a partir de families constitutes de cubes de traces sismiques obtenus 

a Tissue d*une campagne sismique 3D, chacun correspondant k un mSme angle d' incidence. 

6) Methode selon la revendication 1 ou 2 ou 3, caract6risee en ce que Ton forme des 
variables synth6tiques a partir de families constitutes de cubes d'attributs sismiques 
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obtenus par «utem=nt des amplimdes de m>ecs sismiques ob«nu« i Tissue d'une 

campagne sismique 3D. 

7) M6thode selon I'une des revendications prec6dentes, caract6ris6e en ce que I'on 
forme des variables synthetiques a partir de famiUes constituees de cubes de traces 
sismiques obtenus par des explorations sismiques successives de la zone. 

8) Methode selon I'une des revendications pr6c6dentes, caract^iisee en ce que I'on 
fonne des variables synth6tiques a partir de families constitu6es de cubes d'attributs 
sismiques formfis par traitement des amplitudes de traces sismiques obtenues par des 
explorations sismiques successives de la zone. 

9) Mdthode selon la revendication 1 ou 2 ou 3, caract6ris6e en ce que I'on forme des 
variables synth6tiques ^ partir de families constitutes de mesures d'attributs diagraphiques 
de nature physique differente. 

10) Methode selon la revendication 9. caract6ris6e en ce que i'on forme des variables 
synth6tiques h partir de families de mesures de r6sistivite, 

11) Methode selon la revendication 9, caract6ris6e en ce que Ton forme des variables 
synth6tiques k partir de families de mesures lites k la radioactivity des terrains. 

12) M6thode selon la revendication 9, caract6ris6e en ce que Von forme des variables 
synthetiques k partir de families de mesures liees & la porosit6 des terrains. 



# 





FIG.5B 

7*;. so. "ts. i nn- ;^ 



FIG.5C 




250. 2S0. 



• 200., 200. 



ISO. JSD. 



75. 100. SZS. 




soo. 



ISO. El 



too. 100 



50. 50 



15. 100. 12S* 




F1G.6 



FIG.7 




<^00 p 4300 



0 2200 



250 



5/5 



line 



FIG.8 



1^ Class 1 : sand turbidites 
H Class 2 : mud turbidites 
O Class 3 : hemipelagite 



250 



200 



500 meters 



© 150 
100 



50 



J ' I 



r=i: 



T"T 







-4- 


» 1 

— t— 4 — ! 









10 50 90 130 



FIG.9A 

Line 



FIG.9B 

Line 





NATtOMAtDC 

LA ruemm 

UfDOSTnlCLLC 

DEPARTEMEffT DES BREVETS 

-26 bis. rue de Saint P^tersbourg 

*:75800 Paris Cedex 08 

:Tigl6phone : 33 (1) 53 04 53 04T£l«copIe : 



Vos references pour ce dossier (facuUaiif) 


4966/BNI/JG 


D'ENREGISTREMENT NATIONAL 




TITRE DE L'lNVENTION (200 caracteres ou espaces maximum) 

METHODE POUR COMPACTER ET FILTRER DES EVENEMENTS GEOPHYSIQUES DANS LE BUT D'EN 
EXTRAIRE DES INFORMATIONS SUR LA NATURE DU SOUS-SOL 


LE(S) DEMANDEUR(S) : 




INSTITUT FRANCAIS DU PETROLE 

1. ' * * 

i- :. . 




DESIGNE(NT) EN TANT QU'INVENTEUR(S) : 


Q Nom 




Prenoms 


OlMer 


Adresse 


Rue 


1 , aie de la Convention 




Code postal et vllle 


|7i5i0i1 rSI PARIS 


Societe d'appartenance (facullatif) 




B Nom 


FOURNIER 


Pr6noms 


Frederlque 


Adresse 


Rue 


12bis, rue Garnet 


*»',*• . 


Code postal et ville 


19 i5 i2 i4 lO 1 CORMEILLES EM PARIf5Lq 


f - Societe d*appartenance (facultatif) 




; E] Nom 


ROYER 


Pr§noms 


Jean-Jacques 


Adresse 


Rue 


19, rue Beauregard 




Code postal et ville 


15 i4iOiOiOI NANCY 


Society d'appartenance (facultalif) 




S'tl y a plus de trojs inventeurs, uttiisez plusieurs formulaires* Indiquez en haut a droite le de la page suivi du nombre de pages. 


DATE ET SIGNATURE(S) 
DU (DES) DEIVIANDEUR(S} 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 




Alfred ELMALEH. 





l^lPi rf 78-17 du 6 Janvier 1978 relative ^ I'informatique. aux fichlers et aux libert^s s'appllque aux r^ponses faites d ce formulalre. 
Elle garantit un droit d*acc^s et de rectification pour les donn^es vous concernant aupr^s de IMNPI. 

I'. ; 



^ BREVET D'lNA^TION 



33 a) 42 94 86 54 



CERTIFICAT D'UTILITE n- 11235 03 

Code de la propri6t6 intellectueite - Livre VI 



DESIGNATION DMNVENTEUR(S) Page NM . . / 1 . . 

(A foumir dans le cas oii les demandeurs et 
les inventeurs ne sent pas les memes personnes) 

Get imprime est ^ remplir lisiblement d Tencre noire oa 113 @ w/ z/ceoi 



